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Abstract We compute the joint probability density function (jpdf) Pn(M,tm) of the maximum M 
and its position tm for N non-intersecting Brownian excursions, on the unit time interval, in the large N 
limit. For N — >• oo, this jpdf is peaked around M = \/2N and tm — 1/2, while the typical fluctuations 
behave for large N like M — y/2N ex sN^ 1 / 6 and tm — 1/2 oc wN~ x / 3 where s and w are correlated 
random variables. One obtains an explicit expression of the limiting jpdf P(s,w) in terms of the 
Tracy- Widom distribution for the Gaussian Orthogonal Ensemble (GOE) of Random Matrix Theory 
and a psi-function for the Hastings-McLeod solution to the Painleve II equation. Our result yields, up 
to a rescaling of the random variables s and w, an expression for the jpdf of the maximum and its 
position for the Airy2 process minus a parabola. This latter describes the fluctuations in many different 
physical systems belonging to the Kardar-Parisi-Zhang (KPZ) universality class in 1 + 1 dimensions. 
In particular, the marginal probability density function (pdf) P(w) yields, up to a model dependent 
length scale, the distribution of the endpoint of the directed polymer in a random medium with one free 
end, at zero temperature. In the large w limit one shows the asymptotic behavior log P(w) ~ —w 3 /12. 
> 

in 



1 Introduction 



Despite several decades of research, there exist very few exact results for finite dimensional disordered 
systems, beyond mean-field or phenomenological arguments. One important model is the directed 
polymer in a random medium (DPRM) which has been an active area of research in statistical physics 
for the past three decades [TJ. Apart from the fact that it is a simple toy model of disordered systems, 
this problem has important links to a wide variety of other problems in physics, such as interface 
fluctuations and pinning [2J|3], growing interface model of the Kardar-Parisi-Zhang (KPZ) variety [I], 
Burger's turbulence [5], spin glasses [6] or in biological sequence matching problems (7]. The DPRM 
is also a minimal model, yet non trivial, of great interest in the context of disordered elastic systems 
which have found many experimental realizations ranging from domain walls in random ferromagnets 
[5] to wetting on rough substrates [§] etc. 

For concreteness we consider a directed polymer on a square lattice, as depicted in Fig. [1] On each 
site with coordinate (i,j) there is a random energy e^ drawn independently from site to site, which is 
a quenched random variable. Here we consider the ensemble of directed walks of length T (as depicted 
in Fig. [IJ starting from the origin, making step diagonaly left or right, and ending at any point on 
the line depicted on Fig. Q] as the x-axis: this is the DPRM with one free end, or in point to line 
geometry. The total energy E(W) for any directed walk W is just the sum of site energies along the 
walk E(yV) = y~]u ,")cw e v ana - one ^ s interested in the path having minimum energy, that is the ground 
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Fig. 1 Directed polymer of length T on a square lattice. On each bond, there is a random variable dj. 



state or the optimal path. There are two natural quantities describing this optimal path: its energy, 
E opt (T) and the position of the endpoint X(T) [2113] (see Fig. [1]). Earlier studies of this kind of model 
focused on the roughness of this ground state configuration, which is a geometrical characteristic, and 
it was found that 0[3] 



X 2 (T) ex T 4/3 , T > 1 



(1) 



where — stands for an average over the disorder, that is over the random variables e,j. 

Similarly the two first moments of the energy of the ground state E opt (T) were found to scale [2J 
[3J like 



E opt (T)~aT, E 2 pt (T) - E opt (T) ex T 2 / 3 , T » 1 . (2) 

These behaviors ([TJ [2]) were later proved rigorously for one specific model of (discrete) directed poly- 
mer [ 1 OJ . Quite remarkably, more recent studies were able to obtain the full distribution of the energy 
E op t{T) which takes the form, for T > 1 



E opt (T) = aT + b X T 



1/3 



(3) 



where a, b are non-universal constants, depending on the precise model of DPRM under consideration, 
while x is a random variable such that Pr(x < s) = •T'i(s) where ^(s) is the Tracy- Widom distribution 
which describes the typical fluctuations of the largest eigenvalue of the Gaussian Orthogonal Ensemble 
(GOE) [11] , This was shown, rather indirectly, in Ref. [10] using a relation between certain observables 
of the DPRM in point to point geometry and associated ones in point to line geometry [T2] for which 
the limiting distribution had been evaluated [13 . In the context of stochastic growth models, namely 
for the Polynuclear Growth Model in flat geometry, this was shown in Ref. [13]. This was then also 
shown, later, by a direct computation using a relation with non-intersecting Brownian motions [15, 
[T5] and more recently in Ref. [17J using probabilistic tools. Note that in the case where the end of 
the polymer is fixed (in point to point geometry), the same scaling J3]) holds, with different constants 
a',b', and the distribution of \ is a l so different. It was indeed shown [TUl[PH[T51ll9) that in this case 
the random variable \ is distributed according to J-2, the Tracy- Widom distribution corresponding to 
the Gaussian Unitary Ensemble (GUE) [20J, establishing also a connection with the seminal paper of 
Baik, Deift and Johansson [21] on the longest increasing subsequence of a random permutation |22) . 
Recent approaches have also studied continuum models of directed polymers both in point to point 
[231I24U25U26P27U281I29U301I3T] and in point to line geometries [321133] . 

For the directed polymer in point to line geometry, i.e. with one free end as in Fig.[TJ it was further 
shown [TU] that the fluctuations of E opt — E opt and X are described, up to a non-universal rescaling, 
by the statistical properties of the following process Y(u) 



Y(u)=A 2 (u)-u 2 
where A2(u) is the Airy2 process [T§]. Indeed, one has 

T- ' 



lim 

T->oo eg 



(E opt (T) - E opt (T)) = maxy(u) = m 

u 6 R 



M 



(5) 



7V = 4 




Fig. 2 Cartoon of N non-intersecting Brownian excursions on the unit time interval. Here we study the jpdf 
Pn(M, tm) in the large TV limit. 



while 



lim 

T-s-oo 



T 



X — argmaxy(u) = t 



(6) 



which is simply the position at which this process Y(u) (|3]) reaches its maximum. In Eqs. (J3J [5J) ? the 
constant eo and £ are non-universal amplitudes, which depend on the microscopic details of the model 
under consideration. In other words, these relations (JSJ [5]) mean that the joint probability density 
function (jpdf) P(m,t) of the rescaled energy and the rescaled position of the endpoint are given 
by the jpdf of the maximum of the Airy2 process minus a parabola Q and its position. In a recent 
paper |17] , Moreno Flores, Quastel and Remenik have computed this jpdf using rigorous probabilistic 
tools concerning determinantal point processes. They obtained an expression [see Eq. (|127[) below] in 
terms of the product of T\ and a double integral involving in particular the integral operator B m [see 
Eq. (IM|) below], which also enters the expression of T\ (|3"3")l as a Fredholm determinant |34j . 

Here we follow a completely different approach and use the fact that Y(u) also describes the 
fluctuations of the top path of N non- intersecting Brownian excursions in the large N limit [H21l35ll3l)] . 
We remind that a Brownian excursion is a Brownian motion starting and ending at the same point, 
and conditioned to stay positive in-between. We obtain an explicit expression for the jpdf, different 
from the one obtained in Ref. [17 as the product of JF\ and an integral involving a function which 
can be expressed in terms of the so-called psi-function [3"8~1[3"5] for the Hastings-McLeod solution of the 
Painleve II equation [see Eq. (fTTj) below]. Following previous works [15,41,42,43], we thus consider the 
problem of N non-colliding Brownian motions conditioned to stay positive on the unit time interval (see 
Fig. HI) 



< Xi{t) < ■ ■ ■ < x n (t) , V t <= [0, 1] 



(7) 



We focus on special configurations called "watermelons with a wall" where the walkers start, at time 
r = and end, at time r = 1, at the origin (see Refs. [35 , 44J for a rigorous definition of this model). 
We focus on the maximal height M and the time at which this maximum is reached (see Fig. [5]): 



M = max xn(t) , xn(tm) = M . 

0<T<1 



(8) 



The distribution of M was computed for N = 2 in Ref. [33] and then for any value of N in Ref. [31,45, 
|4"5] , Remarkably, it was shown [T5] that the cumulative distribution of M, Fn(M), coincides, up to a 
prefactor, with the partition function of Yang-Mills theory on the sphere, with the gauge group Sp(2/V) 
]47,4S] (see also Ref. [31] for the relation between non- intersecting Brownian motions and gauge field 
theories). Hence Fn(M) exhibits a third order phase transition [501l51ll5"2"] for M = \/2JV, separating 
the left tail of the distribution (which corresponds to the regime of strong coupling in the Yang-Mills 
theory) from the right tail (which corresponds to the regime of weak coupling in the Yang-Mills theory) . 



In the critical regime, for M close to y/2N, which is described by a double scaling limit [47,53], it was 
shown that Fn(M), appropriately shifted and scaled, converges to T\ [TS] [see also Eq. (|3"Tjl below]. A 
rigorous proof of this convergence was recently given in Ref. [16] using Riemann Hilbert techniques. 

The jpdf of M and tm ©, Pn(M, tm) in Eq. ([Hi]) , which is the starting point of our analysis, was 
computed in Ref. [421143] . In this paper we provide a large N analysis of this jpdf. Our main results 
can be summarized as follows. We show that 

lim 2-iN~ip N (V2N + 2- 1 s- sN-i,- + 2-%wN- : 5) = P(s,w) , (9) 

JV— >oo 2 

where the jpdf P(s,w) is given by [J 

4 f 00 

P(s, w) = — J-"i(s) / h(x,w)h(x,—w) dx , (10) 

7T 2 Js 

where -7-1 (s) is the Tracy- Widom distribution for /3 = 1 [see Eq. ([32]) below] and h(s, w) is given by 

h(s,w)= / ($ 2 ((, S )e- w < 2 d(, (11) 

Jo 

in terms of ^(C s) which is one of the two components of the psi-function for the Hastings-McLeod 
solution to the Painleve II equation [3"51l39| . To define this psi-function one considers a Lax pair asso- 
ciated to the Hastings-McLeod solution of the Painleve II equation [35], i.e. the following system of 
linear differential equations for a two-dimensional vector 'P = <^(£, s), 

—& = A&,—& = B&, (12) 

oQ os 

where the 2x2 matrices A = A(£, s) and B = B((, s) are given by 

*<••> = (- 4 c-f-V + 2 / <1 + - + 4cT +2 ■ fl «-)"(-c- f .) • < 13) 

with q = q(s) and r = r(s), where q(s) is the Hastings-McLeod solution of the Painleve II differential 

equation 

q"(s) = 2q 3 ( S )+sq( S ) , q(s) ~ Ai( S ) , (14) 

s—^oo 

where Ai(s) is the Airy function |40) . These matrices A and B constitute a Lax pair associated to 
q(s): this means that the compatibility equation of this system (fT2"|) , i.e. dgd^ — d^dgty, is that q(s) 
satisfies the Painleve-II equation (J14I) and that r(s) — q'(s). The function ^(C: s ) which enters into the 
expression of the jpdf P(s, w) (jlOllip denotes one component of the unique solution (<Pi (£, s), <P 2 ((, s)) 
of the Lax pair (fT2"|) which satisfies the real asymptotics 



*i(c, *) = cos Qc 3 + sA + oic 1 ) , mc s) = - sin Qc 3 + <) + ©(r 1 ) 



(15) 



as C — > ±oo for s real |39j . 

The marginal pdf P(w) is given by integrating the jpdf P{w, s) over s. Although the obtained 
explicit expression for P(w) remains complicated, one can study its asymptotic behavior for w — > oo 
from Eqs. (jlOllip and show that 

1"°° 1 
P(w) = P(s,w)ds , log P(w) = -w^ + oiw 3 ). (16) 

As mentioned above, the results which we obtain here for the extreme statistics of N non-intersecting 
Brownian excursions in the large N limit allow to compute the jpdf P{m, t) of the maximum of the 

1 The factor 2~ 9 ^ 2 in the formula © ensures the normalization of the jpdf P(s,w). 

2 To avoid any confusion we have used a notation for this psi-function which is different from the one used in 
Ref. (1511551175] . Indeed $i and & 2 that we use here actually correspond to & 1 and $ 2 used in Ref. [ToTl39ll73| . 



Airy2 process minus a parabola Y(u) ^ and its position [Ti5, 35, 42, 43] . As explained above, this in 
turn yields important information concerning the directed polymer with one free end (|5I6I) . One expects 
indeed [251131] 

P(m, t) = a(3 P(am , fit) , (17) 

and we show that a = 2 2 / 3 , (3 = 2 4 / 3 . From Eq. (|TB|) . together with ([T7)l . one finds that the marginal 
pdf of the position (O P(4) decays like 

P(i)-exp(-^i 3 j , t»l. (18) 

This result (fT8"|) establishes on firmer grounds a long standing conjecture [1], based on a scaling ar- 
gument and on the study of an approximate model, the so-called "toy-model" [54] [where the Airy2 
process in Y(u) ^ is replaced by a Brownian motion (see also Refs. [55,56])], yielding logP(i) ex — t 3 . 
This behavior is also consistent with numerical studies [42U57B58] . 

As it is well known, directed polymer models as in Fig. Q] can be mapped onto stochastic growth 
processes in the KPZ universality classes in 1 + 1 dimensions [1,59]. It was shown in particular [14,19] 
that the fluctuations of the height field, in such processes, are governed by T<i in curved (or droplet- 
like) geometry, and by T\ in flat geometry. Recently, these theoretical predictions for KPZ interfaces 
have been observed in remarkable experiments on turbulent liquid crystals, both in curved and in flat 
geometries |60,:6I]. The quantities which we focus on here (J3J EJ) have also a clear physical meaning for 
growth processes in the droplet geometry [42"ll43| : the maximum of Y{u) corresponds to the maximal 
height of the droplet, while its position corresponds simply to the position of the maximal height of the 
droplet. Very recently, the distribution of the maximal height of this droplet was measured in the same 
experiment on turbulent liquid crystals and a very nice agreement with the Tracy- Widom distribution 
T\ [see Eq. ([ST]) ] was found [62], as expected from Eqs. (J3J[5j. 

The paper is organized as follows. In section 2, we give an expression of the jpdf of M and tm 
for TV excursions (see Fig. ^ in terms of discrete orthogonal polynomials. In section 3, we analyze the 
large N and large M limit, when M is much larger than its mean value, M ~^> y/2N. In section 4 we 
analyze the system of orthogonal polynomials in the double scaling limit, which allows to compute the 
limiting jpdf P(s, w) (|9"1 [TU]) . before we conclude in section 5. We have left in Appendix A the detailed 
analysis of the tail of the marginal pdf P(w). 



2 The joint probability density function in terms of orthogonal polynomials 

In this section we derive an expression of the jpdf P^(M,tm) in terms of discrete orthogonal poly- 
nomials, which turns out to be very useful to perform the large N analysis. The starting point of our 
analysis is an exact expression, using path integral for free fermions [411163] . for the jpdf Pn{M,tm) 
of M and tm given by [4"2"Il4"5] 



P n (M,t m ) = 9N+lMN N (2N+1)+3 J2 



n,n' N 



JV-1 



(_l)""+»i, n %n' N 2 H nlA N (nl...,n 2 N _ 1 ,n 2 N ) (19) 



i=l 



N-l 



xA N (nl...,n 2 N „ 1 ,n' N 2 )e ^ & "VsSM* 1 -™ Kv 2 +W, 

where Z\jv(Ai, . . . , Ajv) is the N x N Vandermonde determinant, and where we use the notations 
n = (rii, . . . , rtjv) and ^ n = '^2 n= _ 00 - The numerical constant An is given by |43| : 



Ntt 2N + n + 2 



2^-^nf=o 1 r(2+j)r(i+j) ' 



(20) 



We thank J. Krug and J. Rambeau for pointing out these references [571158] . 



which ensures the normalization of Pn{M,tm) |43| . i.e. 



dM dr M P N {M, TM ) = l 
o Jo 



(21) 



To study this multiple sum (|19l) we introduce discrete orthogonal polynomials pk(n) such that |47) 

oo 

2 , Pk{n)pk>{n)e ^^ = 5 k ,k>h k , (22) 

n=— oo 

where Pfc(n)'s are monic polynomials of degree A;: 



and /ifc's are positive constants. In particular one has 

N 



N 



YlmA{nl 



2 „2 



, u 2 , ■ ■ • , n N ) 



det pai-ifnj) = ^Z e (°") T]P 2CT M-i( r 



(23) 



(24) 



o-e5i\ 



where <S/v is the group of permutations of size TV and e(cr) is the signature of the permutation a £ iSjv- 
After some manipulations, using this expansion ([24)) together with the orthogonality condition (j22|) we 
arrive at the following formula for Pn(M,tm)'- 



N 



P N (M,T M ) 



(N-1)\A N 
2 N + 1 M 2N ' I+N + 3 J- J-" 2 "' ! 

3=1 



n** 



V(-l)" +m nTO V fe^Wfe-'H ^^ld-rMl^+rMm 2 ] _ 



h 2 i-l 



If one defines tjvj = | + u one finds Pjy(M, tm) = Pn(M, u) with 



(25) 



P N (M,u) 



(n-i)\a n 

2 N + 1 M 2N2+N+3 



n 



N 



Jl h 2j -i Y, G2k-i(M, u)G 2k -i(M, -u) 



3 = 1 



fc=l 



G 2k -i(M,u)= ]T (-l) n n^-!(r 



n— — oo 



where we have used the standard notation ^fc(n) for the "wave function": 

, i pk(n) -^fL„ 2 
V^fc(n) = e Im 7 



(26) 

(27) 

(28) 



Using the result of Ref. [15], one notices that Pn(M,tm — 1/2 + u) in Eq. ([26]) can be rewritten as 



JV 



P N (M, u) = F N (M) ^3 ^T G 2fc _i(M, u)G 2fe _i(M, -«) , 

JV! 



fc=i 



F N (M) =Pr 



max xn(t) < M 

0<T<1 



r 2Af 2 +Af 



AT 



nS 1 r ^ 2 + ^ r d + •?') 2 " 2+fMM!+ " ii 



n^ 2 



(29) 



(30) 



where Fn(M) is thus the cumulative distribution of the maximal height of N non-intersecting excur- 
sions, on the unit time interval. In Ref. [15] the large N asymptotic of Fm{M) was carried out and it 
was shown that (see Ref. [TB] for a rigorous proof of this result) 



lim F; 

TV-s-oo 



V(V2AT(l + .s/(2 7 / 3 JV 2 / 3 )) = ^(s) , 



(31) 



where ^i(s) is the Tracy- Widom distribution for {3 = 1 (f3"2]l . We recall that J 7 ! admits the following 
explicit expression [11 

r 1 (a) = exp(-±J {(t- s )q 2 (t) + q(t)) dt) , (32) 

in terms of q(t), which is the Hastings-McLeod solution to the Painleve II equation (|14p . Note that 
.Fi(s) can also be expressed as a Fredholm determinant as |34| 

^ 1 (s)=det(/-77 B s 77o), (33) 

where B s is an integral operator with kernel 

B 8 (x,y)=Ai{x + y + s), (34) 

where Ai is the Airy function and IIq is the projector onto the interval [0, +00). 

The main result of this section is that we have reduced the analysis of the large N analysis of 
Pn(M,tm) to the analysis of the sum over k in Eq. (|2"§|) which involves the quantity G2k-i(M,u) in 
Eq. (EZD- 



3 Analysis of G-zk— i(Af, u) for large M: large deviation analysis 

In this section we study in detail the behavior of G2k-i(M, u) in the limit of large M . More precisely, 
we assume here that M ~^> \/2N, and more precisely 



M - V2N = 0(V N) , (35) 

which corresponds to the "right tail" of the distribution of the maximum (l3"Tj) . In this limit where M 
is large, the discrete sum which defines the orthogonal polynomials pk in Eq. (|22p can be replaced by 
an integral and therefore the orthogonal polynomials pk are well approximated by Hermite polynomi- 
als g7]. For M > V2N one thus has 



p k (x)=(^r) H k (^-x)+0(e-^), (36) 



k 

\V2n) " K V\/2Af 
where Hk is the Hermite polynomial of order k [65 



I - 1 
Hk { z ) = k\ Y { l ] n [Zz l - , (37) 



where \x\ is the largest integer not greater than x, while the amplitude hk in (f2"2"| is given by (see 
Ref. ISj) 

h k = V2^k\(^\ +0{e~ 2M2 ). (38) 

From Eqs. (|3"6")l and (|3"8"|) one obtains that tpk in Eq. I|28p is given by 



With this expression (|3^|) . it is then possible to study G2k-i(M,u) in the limit M ^> \/2N and N 
large. We first consider the case u — (for illustration) and then present the analysis for arbitrary u. 
This then allows us to study the jpdf (I2!J|) in the large deviation regime. 



3.1 The case u = 

To analyse G 2k -i(M, u — 0) in Eq. (|27p it is convenient to use the Poisson summation formula to 
obtain 

oo oo 

G 2k -i(M,u = 0) = J2 (-l) n n^2 fe -i(n)= ]T F(n) , (40) 



n=— oo n— — oo 



nw) = ^ y dx ^~ iyx x fe-i (^j . (41) 

Using the asymptotic behavior of ip 2 k-i(%) in Eq. (|39[) one thus obtains that, for large M, G 2k -i {M, u = 
0)~[F(0) + F(1)]/(2tt), yielding 

Ifi 23~ fe f 00 

G 2fe -i(M,u = 0)~ 7/ . __ =M* / dzcos(2Mz)e-* z H 2k ^(V2z) . (42) 

tt' /4 V^(2fc - 1)! Jo 

This integral can be performed by expanding the Hermite polynomial as |65j 

_ fc -! nf(2fc-2p-l) 

tf 2fc _ l( ^) = (2fe - 1)! g ( - 1) %!(2 fc -2 P -l)! z2fc " 2P " 1 ' (43) 

and then integrate over z term by term in Eq. (|4"2"| . After some algebra one finds that G 2k -\(M, u = 0) 
can be rewritten in terms of Hermite polynomials as 

G 2k . 1 (M,u = 0) - -J— - 4====Afi e - M2 (H 2k (V2M) - MV2H 2k _ t {MV2)) . (44) 

7T4 2 fc v /(2A; — 1)! \ / 

As we show below, this analysis can be extended to any finite value of u. 



3.2 The case of finite u 

For finite u the Poisson summation formula (|40[) yields the following asymptotic expression 



G 2k -i(M,u)^(-l) k —— / rfzcos -==z e z z # 2 fc-i W — — r- « 



7T 7 / 4 ^(2A: - 1)! 1 + 2u 7 " \^T+2^J ^ ' \\ 1 + 2u 

(45) 
Expanding again the Hermite polynomial as before (|43p and integrating term by term, one obtains 

, 2 11 / 4 1 Mle"^ /l-2u\ fe 

"■H-^ >v ;.,i ;v (ri) ,46 » 



' 2U H 2k ( M «/ - 2 , „ | - V2MH 2k -i [ M, 



1 + 2u I V 1 - 4u 2 / \ V 1 - 4u 2 

Given this expression (05]), it appears that the sum over k which enters the expression of Pn(M,tm) 
in Eq. (|29p is actually dominated, for large M, by the large values of k, which we now study. 



3.3 Asymptotic limit of G<2k-i{M, u) for large M and large k 

In the limit of large M we are thus interested in the behavior of G<zk-i{M, u) for large k ~ O(N) with 

fixed. Here c < 1, corresponding toAf> y/2N [see Eq. (1551) ]. To study the behavior of G2k-i(M, u) in 
Eq. (|4"u]l in the limit of large M and large k it is convenient to use an integral representation of Hermite 
polynomials (see also Ref. [471151] for a similar calculation with different orthogonal polynomials). A 
convenient one is the one coming from the relation for the exponential generating function [40 , 65 

°° t p 

2^JJ p (x)-7 = exp(2x£-£ 2 ) , (48) 

p=Q P ' 

which allows to write 

*#>-*} &*****' < 49 > 

where the contour must encircle the value t = 0. From this representation (|49p one obtains 



H "[ M \lrk?) ^ 2fc ) ! /^i exp ( 2 VT^ Mt - <2 - clog V ' (50) 

which is quite convenient for a large M analysis. Performing the change of variable t = yM, one gets 



(2k)\ i dy 1_, ,,, 



H 2k [ M,/ T - i -H = ^£f ^-exp(-M^)) , (51) 



<t>{y) = y 2 - 2 J y~^ v ' + clogy ■ (52) 

For c < 1, the contour integral in (fSTj) can be evaluated for large M by a saddle point, which yields 



f7n.. -, fM 7/^ ~ V 1 






ZK ll ' J l j ^/^^.^(i^u^l^J 




., (2fc)! c -MW, 1 / 1 /l-2« V2\ 


(53) 


M 2fe V2^a/ m2 I^"(2/*)I V y * ^ l + 2u c y 
with 


, 1-Vl-cp ,, „. 2-2Vl-cp + cp(l + ln(2p)-21n(l-Vl- 
y y/Tp ' W J 2p 

rh"U,*\ -9 9 Cp <r n 


-cp)) 


(54) 
(55) 


(1-Vl-cp) 2 
in terms of 



p = 1 - Au 2 . (56) 

On the other hand, using Stirling's formula, one obtains, in the large M limit, keeping c = 2k/M 2 
fixed 67D 



^jS ~ (2 TC ) 1 /4VM e -M^ f+S1 „ 2 - 1 „ c ) (5?) 
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Finally, combining Eq. ([53"]) together with Eq. I|57|l one obtains 

M 2 c 3 / 4 



2 5/2 

G 2k -i(M,u) ~ (-l) fe 



7T3/2 (1 - 2 U ) 3 / 2 Vl0"(r)| I ^ 



1 /l-2u \/2 



l + 2u 



x exp 



"^(■-'"W^'-fS 



1 



<Xy* 



l + 2u 
Let us analyse this formula for u = 0, where G2k-i(M,u = 0) takes a simpler form given by 

t 4 



(58) 



Ga fc _i(Af,u = 0)~(-iy 



^v 1 -^, , 



l-c 



1/4 



x exp 



In the limit where c — > 1 . one finds 



M 2 [ yi - ^ + cln (l - \/T^c) 



1 



-cine 



G 2fc _ 1 (M, U = 0) ~ (-1)*-^M 2 (1 - c )V4 e -|M^(i_c)3 



In particular, in the regime where 



(1 - c) = iir 4 / 3 , 



(59) 



(60) 



(61) 



which corresponds precisely to the tail of the double scaling regime to be studied in section 4, one finds 

G 2k -i{M,u = 0) ~ (-l) fe ^Af 5/ V /4 e-^ 3/2 ~ (-l) fe+1 -Af 5/3 Ai'(x) , (62) 

where Ai'(a;) is the derivative of the Airy function [3D] and where the last relation in Eq. (|6"2"|) holds 
only for large x. Note that this result (|6"2"|) . involving Ai'(x) can also be directly obtained from the 
above expression in terms of Hermite polynomials given in Eq. (|44p . To do so, we make use of the 
Plancherel-Rotach formula [6l)ll66| 

cxp(-z 2 /2)H 2fc+m (z) = (4k)^nh k +i^2ky.(2k)-^ ( Ai(t) - -^rAi'(i) + 0(k^)] , (63) 

V (2k) 3 J 

valid for large k, where we have set 



z = (4fc) 1 / 2 +2- 2 / 3 fc- 1 / 6 i 



(64) 



Indeed, if one applies this Plancherel-Rotach formula (JS3J) to Eq. (gj) with z -> y/2M = (Ak) 1/2 + 
2-2/3^,-1/6^ j gee gq g ^ pyj IDT]) ], one obtains precisely the above asymptotic behavior (|62p in terms 
ofAi'(x). 

It is instructive, and useful for the forthcoming double scaling analysis, to apply this Plancherel- 
Rotach formula (J63|) to G2k-i(M, u) for finite u, as given in Eq. (|46|) . where we also rescale u accord- 
ing to 



v M^ 3 



(65) 



One can then apply this formula $$ to flU with z -» M^2/(l - u 2 ) = (4fc) 1 / 2 + (a;+4u 2 )2~ 2 / 3 fc- 1 / 6 , 
which yields 



G 2fe _i(M,u) ~ (-1) 



fc+i 



^+ 2 ^ (2wAi(4 w 2 + x) + Ai'(4w 2 + x)) 



(66) 



Note that the same function (|6"r!)) . albeit with different arguments, enters the expression of the jpdf 
given in Ref. [17] [see Eqs. (I127p . (11281) below]. Here this function arises naturally as the Plancherel- 
Rotach asymptotic of our large N, large deviation, regime. This fact will allow us to identify the 
coefficients a and /3 in Eq. (1171) . 



3.4 Large deviation regime of the joint probability density function for M ^> V2iV 
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From the asymptotic behavior obtained in the previous section (|58|) it is now possible to obtain the 
behavior of the jpdf P N (M 7 t m ) = Pn(M, u) with t m = 1/2 + u. From Eq. 1^5]). using that F N (M) ~ 1 
when M ^ V2N (see Ref. [15J for a more refined analysis of this right tail) one has 

2 N 

P n (M,t m ) ~ ^^G 2k -i{M,u)G 2k _ 1 {M,-u) ( 67 ) 

k=l 



N 
lvi ^— \ ivi 1 ill — ci ., i zk; 

— L a _ 4u 2)3/ 2 U"(„*)| C l/2 eXP L~ M ^ (C ' U) ] ' °= W 2 (68) 



2 4 M<A M 4 1 2(1 - c) r 2 .-, 2fc 



, (l-4u 2 ) 3 /2 |^"(j/*)| C l/2 



= 2 ^ /1 CP -cln(cp) + 2cln(l- */l - cp) , /? = 1 - 4u 2 . (70) 

P 

One can easily check that for any value of u 6 [—1/2,1/2], ip(c,u) is a decreasing function of c. 
Therefore the sum over k in Eq. ([^7)1 is dominated, for large M, by k = N. This yields, in this regime 
(in the sense of logarithmic equivalent, with 2N/M 2 fixed) 

P N (M, u) ~ exp [-M 2 ip(2N/M 2 , u)] . (71) 

In particular, when M 2 — » 27V, one has 

P N (M, u) ~ exp (-2/V^(l, u)) , (72) 

¥>(1, «) = TT^TI - ln (1 - 4 « 2 ) + 2 In (1 - 2|u|) = f \u\ 3 + 0(\u\ 5 ) . (73) 

This cubic behavior for small u implies 

P N (M ~y/2N,u~vM- l / 3 )~exp(~\v\ 3 ) , (74) 

which exhibits an interesting non-Gaussian behavior. We will recover this cubic behavior (J74I) below, 
in the double scaling. 



4 Double scaling regime 

The above results concern large deviations: they describe the fluctuations for which M is much bigger 
than the mean value, when M ^> y/2N. However, it was shown in Ref. [T5] that the typical fluctuations 
of M behave like M — \^2N — 0{N^ X / 6 ). In the language of Random Matrix Theory, this corresponds 
to a double scaling regime. We thus set M — \/2N = sTV -1 / 6 where s is fixed and we are led to study 
the quantity G 2 fc-i(Af, u) where both k ~ N and M 2 ~ 2N. To study this regime, we need to study 
the discrete orthogonal polynomials p^s (|22p beyond the approximation where p/t's are estimated 
by Hermite polynomials (|36p . To perform this analysis, we study in detail the three terms recursion 
relation satisfied by these orthogonal polynomials j^T]. This recursion relation for p^'s yields differential 
recursion relations for G2k-i(M, u), which we derive in the first subsection. In the second subsection, 
we remind some useful results obtained by Gross and Matytsin |47j concerning the analysis of this 
three terms recurrence in the double scaling limit. In the third subsection we obtain a differential 
equation and a partial differential equation satisfied by C?2fe— \{M,u) in the double scaling limit and 
obtain the limiting form P(s, w) of the jpdf Ppj(M, tm) for large N. In the fourth subsection, we show 
how to solve the aforementioned differential equations. In the last subsection, we present an asymptotic 
analysis of the marginal pdf P(w). 



12 



4.1 Differential Recursion Relations 



In this subsection, we derive two distinct differential recursion relations satisfied by G 2 fc_i(M, u) which 
will be very useful to perform the large N analysis of the expression above (|26[) . For this purpose we 
remind that the discrete orthogonal polynomials pk 's defined by the orthogonality relation ((22)) satisfy 
the recursion relation 



xpk(x) = p k +i(x) + R k Pk-i(x) , R k = , (75) 



such that t/Vs in Eq. (128H satisfy 



xipk(x) = 'f k +iipk+i(x) +jk*Pk-i{x) , 7fc = \f~Rk~ = \ t • (76) 

From this recursion relation (|76p one deduces easily the following differential recursion relation 

d ir 2 

-Q-Gzk-x = ~~^TH [72fe72fc+iG 2 fe+i + {"flk + llk-i)G2k-i + 72fc-272/£-iG 2 /c-3] j (77) 

where, for clarity, we used the notation Gik-\ = G2k—i(M,u). 

One can also derive a differential recursion relation for G2k-i(M,u) when M is varied. By differ- 
entiating the orthogonality relation (|22[) with respect to M, one obtains the identity (see also Ref. [69J 
for a similar relation for more general matrix models): 

M 3 d ( I , 1 

- Y~5~dM^ k = X lklk ~ 1 ^ k ~ 2 ~ 4^+27/0+1^+2 , (78) 

from which one obtains straightforwardly: 



M 3 d 

2^dM 2k - X 



1 U\ U , 2 2 . 

4 _ 2 J 72fe-i72fe-2C J2 fc-3 - 2 (72fc + 72fe-iJ G 2 k-\ 



1m, 

T + 7j J 72fe72fc+l ( J"2/c+l 



(79) 



These two differential recursion relations (|77l 179")) will be analyzed below in the double scaling limit 
where they will lead, respectively, to a partial differential equation and an ordinary differential equation 
for G 2 k-i- 



4.2 Three terms recursion relation in the double scaling limit 

In this regime, it was shown by Gross and Matytsin [47] and later, rigorously, by Liechty [16], that the 
coefficients Rk in Eq. (|75p take the following scaling form: 

R2k = ^ - M 10/3 /i(x 2fc ) + M^f+{x 2k ) + 0(M 2 ) , x 2k = M 4 / 3 (l - |^) , (80) 

R2k+i = -~2-+ M 10/3 fi(x2k+i) + M 8/3 f 2 -(x2k+i) + 0(M 2 ) , x 2k+ i = Af 4 / 3 M - -j^j ,(81) 

where the function f\ satisfies a Painleve II equation (PII), corresponding to a = 0, 



7T 



2 



o5/3 

f 3 /^ f.w ~ a,-/"? 2 / 3 



/{'(a:) = 4a:/i(x) + —ff(x) , A(x) T M{2^ A x) , (82) 

Z £^-00 7T 
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where Ai(y) is the standard Airy function. Note the minus sign in the asymptotic behavior wich is 
missing in Ref. |47| . It can be expressed in terms of the Hastings-McLeod solution of PII given in 
Eq. Hm as 



2 5/3 



h{x) = r q(2^ 6 x). 

One can also show [57] that the functions f 2 , f% m Eq. ([50)1 satisfy 

f+(x) + f 2 {x) = -JL x + £f?(x) . 



(83) 



(84) 



We also notice that f 2 = f 2 while this identity is not really needed here (in fact one can check that 
this is always the combination f 2 + f% which enters into the calculations [70]). 



4.3 Recursion relations for (j?2k-i(M, u) in the double scaling limit 

From the large deviation analysis performed above and in particular from the behavior obtained in 
Eqs. (IHll); (|66p in terms of the variables x in ([ST]) and v in (|65p . one also expects that G 2 k-i(M,u) will 
be, in the double scaling limit, a function of the variables 



r,,= U' :i (l--^) . r -uM^'- 



(85) 



On the other hand, due to the term (— l) fc+1 in (|55|) . one expects that that G 2 k~i(M, u) behaves 
differently for k odd or k even, also in the double scaling limit [47] . Therefore, guided by this analysis, 
and by the behavior of Rk in Eq. (|50|) , one assumes the following ansatz 

G 2k -i(M,u) = M 5 / 3 (g+(x 2k ,v) + M- 2 / 3 g+(x 2k ,v) + M- 4 / 3 g+(x 2k ,v) +0(M- 2 )^j , k even(86) 

G 2k -i(M,u) = M 5 ' 3 ^(x2k,v)+M-y 3 g2 -(x2k,v)+M-^ 3 g^(x2k,v)+0(M- 2 fj ,fcodd,(87) 

where gi , g 2 , g 3 are some functions which remain to be determined. From Eq. f)80[) one has 

72fc72fc+i - y/RzkRzk+i = ^T- M " /3 f^/i (a*fc) + ^/{M + ^safc") + 0(A/ 2 ) 

7| fe +72 2 fe-i=^2fe+i?2fc-i = 2. 



-^ + m 8 / 3 (y A 2 (x 2fc ) + /{(a**) - A* afc ) + 0(M 2 



(89) 



.5.i The case u = 



Substituting this ansatz (|5rJl) into Eq. ([79"]) . for u = 0, and performing the same analysis as the one 
done in Ref. [37] one obtains 



gf(x, 0) = -g 1 (x, 0) = gi(x, 0) , 
where <?i(x, 0) is the unique solution of the third order linear differential equation 



(90) 



—y'"(x)+y'(x) 



ui(x) = 



1 



!*/ 



/!(*) 



+ ^A 2 (*) + 



- 


!/(*) 


a; 




7T 2 _ 


) 



^+4ui(a:) 



whose large x behavior has to match the large deviation regime, where M 
Eq. 



, (91) 

(92) 

V / 2 _ /V from above [see 



in (x, 0) 



-Ai'(x) , a? — > +oo . 



(93) 
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As a check, if one uses the large x behavior of u\(x) in Eq. (1921) . U\(x) ~ a;/(47r 2 ) one finds that the 
above third order differential equation reads (|9ip . for large x 



y'"(x) +xy'(x) - 2y(x) = 



(94) 



which admits y(x) = Ai'(x) as one of its three independent solutions. It is convenient to write this 
differential equation (|9ip in terms of the Hastings-McLeod solution q(s), which is related to fi(x) 
through the above relation (|83[) . In particular one finds straightforwardly 



u ^ x ) = w^ iq2{s) - q/{s) + i hs = 2V3x ■ 

From Eq. (J9"Tj) one has also 

9l (x,v) = f(s = 2 2 ^x,w^2 7 / 3 v), 
where f(s,w — 0) is the solution of 

V'(«) - V'(s) [6 {q 2 - q') + a] - y(s) [s {q 2 - q')' + 2 

which behaves, for large s, like 

f{s,w = 0) ~ --Ai'(2- 2 / 3 s) , s-> oo . 



= 



(95) 
(96) 
(97) 

(98) 



4.3.2 The case u ^ 

For u 7^ one obtains that <?i(a:, i>), with u = tiM -2 ' 3 is given by the solution of the following 
differential equation where v plays the role of a simple parameter 



4 16 

^V*(x) - v-^Vvfr) + Vvfr) 

-Vv{x) 

whose large x behavior is given by Eq. (|66p : 



2x 
3^ 



8ui(x) 



^ An 

+ iu ' 1 (x)--(7T 2 f 2 (x)+2f[(x)) 

&-K Z 3 



= 



gi (x, v) ~ -- e ^+ 2 ^ (2tAi(4v 2 + x) + Ai'(4t> 2 + x)) . 

We now perform the change of variable 

9l (x,v) = f(s = 2 2 / 3 x,w^2 7 ^v), 
and one then obtains, from Eq. (|I00[) . an equation satisfied by f(s,w) 

\yZ - 2wyl - y' w [6(q 2 - q') + s] - y w [3(q 2 - q')' + 2 - 2w(q 2 - q')] = , 



(99) 
(100) 

(101) 

(102) 
(103) 



where y w = y w (s) and where f(s,w) is given by the unique solution of Eq. (|103l) which, for large s, 
behaves as 



f(s,w) 



. 2 U/3 j_ 

e 24 



1 
2"27 3J 



^Ai( W 2 /2 8 / 3 + s/2 2 / 3 ) + ' Ai'(u; 2 /2 8 / 3 + s/2 2 / 3 ) 



(104) 



On the other hand, from the recursion relation (|77l) . one obtains, in the double scaling limit, that 
f(s,w) satisfies the following equation 



i Hs - w) = 



£-<' s -''» 



f{s,w) , 



(105) 



which has actually the structure of a supersymmetric Schrodinger equation. This in turns means that 
/(s, w) satisfies a Fokker-Planck equation. The connection between the Airy2 process and this diffusion 
process deserves certainly to be explored. This equation f|105[) . together with the third order differential 
equation satisfied by f(s,w — 0) (|103|) determines uniquely the function f(s,w). 
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4-.S.3 Limiting form of the jpdf in the large N limit 
We recall that the jpdf Pjv(M, tm) is given by 

2 N 1 

P N (M,T M ) = F N (M)^ ^G 2fe _i(M,u)G'2 fe _i(M,-u) , r M = ^ +^ . (106) 

fc=i 

using the scaling form found above (|86p . one finds 

P N (M,r M ) = F N {M)^-M 1 ' 3 (107) 

f> (^=^, (tm - 1/2)M 2 / 3 ) * (^*, -(r M - l/2)^/») + <D(M-y*) . 



x 
fc=i 



If we set 

M = V2N + 2- 11/6 N- 1/6 s , i.e. M 2 = 2N + 2" 1/3 7V 1/3 s + 0{N~ 1 / 3 ) , 

T M -- = 2~ 7 / 3 wM-y 3 = 2-^wN- 1 / 3 , (108) 

one finds that for large N, using also the result above (f3"Tj) , obtained in Ref. |15| 

P N (V2N + 2- 11/6 N- 1/6 s, 2- 8/3 N- 1/3 w) (109) 



•2 



N 



where we have simply made the change of variable j = N — k. In the large N limit the discrete sum 
over j converges to an integral where, thanks to the behavior of f(x, w) for large x the upper bound 
of the integral can be safely taken to infinity. Finally one obtains the large N limit as 

lim 2-iN-ip N (V2N + 2- 1 s- sN-^,- + 2-^wN-i) = P(s,w) , (110) 

JV-s-oo 2 

where P(s, w) is given by 

2 roa 

P ( s i w ) = 02073 -^W / f(x,w)f(x,-w)dx > (HI) 

where /(x, to) is the unique solution of the third order differential equation in Eq. (|103p which has 
the asymptotic behavior given in Eq. (|104p for large a;. It also satisfies the Schrodinger-like equation 
in Eq. (|105p . In the following we show how to solve this equation f|103[) in terms of a psi- function 
associated to the Painleve-II equation. This second equation (I105p turns out to be very useful to guess 
the form of the solution. 



4.4 Solution of equations (|97[) and (|105p in terms of psi-function associated to Painleve-II 

Such differential equations (|103l I105|) where the coefficients involve Painleve II transcendents have 
appeared several times in the literature on related subjects (see for instance [T3], [71], [75] ). In these 
cases it was useful to introduce Lax pairs associated to the Painleve II equation. On the other hand, 
in Ref. [T^ it was shown that the Christoffel-Darboux kernel associated to the orthogonal polynomials 
studied here (j2"2"|) is described, in the double scaling limit, in terms of psi-functions related to a Lax 
pair for the Painleve II equation. This gives a hint that this Lax pair might be relevant to analyze these 
differential equations (|1031 1105[) . We introduce briefly this Lax pair and refer the reader to Ref. [TBI 
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[75] for more detail. Following Ref. p&L one considers the linear differential equations for a 2- vector 
«P = «P(C,s), 

— V = A$, — & = B&, (112) 

oQ os 

where the 2x2 matrices A = A(£, s) and B = B((, s) are given by 

A(C,)-( ^ q 4C 2 + S + 2 9 2 + 2r\ 

and 



B(Cs)=( K q c J q j . (114) 

These matrices A and 5 constitute a Lax pair associated to the Hastings-McLeod solution of the 
Painleve II equation: this means that the compatibility equation of this system (|112[) . i.e. d s dc& = 
d^dsty, is that q satisfies the Painleve-II equation and that r(s) — q'(s). We denote by 



*i(C,«) 

*a(C,») 

the unique solution of the Lax pair (J112I) which satisfies the real asymptotics 



(115) 



*i(C> *) = cos (^- + sCJ + OiC 1 ) , *a(C, «) = - ^ (^- + <) + ^(C 1 ) , (116) 

as £ — >■ ±oo for s real. There is such a solution (see [39]) and it further satisfies the properties that 
^i(C, s) and ^ 2 (C> s ) are rea l f° r re &l C an( l s an d 

* X (-C, s) - <Pi(£ s) , *a(-C, s) = -* 2 (C, s) . (117) 

One can also check from (|112[) that, for s — > oo and £ real, one has 

<Pi(£ s) - cos f ^- + <") , <£ 2 (C, *) ~ - sin (^- + sA . (118) 

Note that the psi-functions initially studied by Flaschka and Newell in Ref. [35] are <P± + i^ 2 a nd 
<Pi — i$>2- From Eqs. (J112U114]) . one finds that ^ 2 (£, s) satisfies the interesting equation 

9> 2 (C, a) - (q 2 - <7> 2 (C, «) = -C 2 ^ 2 (C, *) ■ (H9) 

This relation (|119p , in view of the above equation (|105[) suggests to look for a solution of (|103p under 
the form 

/•OO 

f(s,w)= / C (C)<P 2 (C,s)e-^ 2 dC- (120) 



To find this function c(£) it is sufficient to look at the large s behavior of f(s,w = 0). Indeed we 
obtained previously that 

/(s,w = 0)~--Ai'(2- 2 / 3 s). (121) 

On the other hand from Eq. (|120[) and Eq. (|118[) one has that 

f{a,w = Q)~-J c(0sm(^ + S (]d(. (122) 
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Now using the integral representation of the derivative of the Airy function [40] 



one recognizes that 



i r°° /i 3 

Ai'(z) = / isin I — + zt ) dt , (123) 



oi3/3 
c(C) = -C- (124) 



One thus finds that f(s,w) can be expressed in terms of ^(C s ) as 

2 13 / 3 f 00 ,2 

f(s,w) = {<Z> 2 ((,s)e- wi d(. (125) 

71 Jo 

We can then check, using the equations in (|112[) that this expression given in Eq. (]125[) is indeed a 
solution of the third order differential equation given in Eq. (|103[) and behaves asymptotically as in 
Eq. (|104j) . This expression (|125p together with the Eq. (|llip yields our main result announced in the 
introduction in Eqs. ([10]). (|TT| . 

Looking at the expression (I125p one might be worried about the convergence of this integral for 
w < 0. However, when £ — > 00, ^(C s ) 1S a highly oscillating function, ^(Ci s ) ~ — sin(|£ 3 + sQ, and 
it is then possible to give a meaning to this integral by using a standard regularization scheme, where 
one multiplies the integrand by exp (— e|£| 3 ), with e > and then take the limit e — >• + . One can 
check explicitly this procedure to obtain that the large s behavior of f(s,w) as defined in Eq. (]125[) 
yields the expected behavior in Eq. (|104l) . In that case one indeed replaces ^(Ci s ) hy — sin(|C 3 + s(), 
[which is its large s behavior (I118p ] and the integral over £ in Eq. (|125p can indeed performed, also for 
w < 0, using the aforementioned regularization. 

4.5 Asymptotic analysis 

Here we first give the asymptotic analysis of the joint pdf for large s. As we will see, this provides a 
way to compare our results with the results of Moreno Flores, Quastel and Remenik given in Ref. |17| 
and to obtain the amplitudes a and /3 in Eq. (|17p . Then we obtain the asymptotic behavior of the 
marginal pdf P(w) for large w. 

4.5.1 Joint probability density function for s — > 00 

We first begin by analyzing the limit s — > 00 of the jpdf P(s,w) in Eq. (jllip . In this limit, one can 
obtain an estimate of the integral over x in Eq. (jllip by simply replacing h(x,w) by its asymptotic 
behavior given in Eq. (|104[) as 3; > s > 1. Besides, for s — > 00, one has to leading order IFi(s) ~ 1. 
Therefore one has, for s — > 00 



P(s, w) ~ / dz 

/2- 2 /3 s 



oc 



Ai ' 2 fi+-)-!^ Ai2 (x+- 



(126) 



On the other hand, in Ref. [17], the authors obtained the following expression for the jpdf P(m, t) 
of the maximum $5§ and its position © of the process Y(u) ([4]) [see Eqs. (1.3), (1.4) therein]: 

P(m,t) = 2 1 / 3 J' 1 (2 2 / 3 TO) / dx dy$_ t:m {2 1 / 3 x)p 2 2 /3m (x 7 y)<P t M2 1/3 y), (127) 

Jo Jo 

where 

$ t ,m(x) = 2e xt [tAi{t 2 +m + x)+ Ai'(t 2 + m + x)] , (128) 

and 

p m (x, |/) = (J- noBrnllo)- 1 ^, y) , (129) 
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where B m is the integral operator defined in Eq. (j3~3j) and IIq is the projector on [0, oo). In the large 
m limit, one can use the estimate, p 2 2/3 m (x,y) ~ 5(x — y) in Eq. (|127[) . so that one gets the large m 
estimate of P(m, t) as 

/>oo 

P(m, t) ~ 4 / dz [Ai' 2 (f + z) - i 2 Ai 2 (i 2 + z)] . (130) 

By comparing these two asymptotic formulas (|126l 11301) , one finds that a and /3 in Eq. I|17|) are given 
by a = 2 2 / 3 and /3 = 2 4 / 3 such that 

P(m, t) = 4F(2 2 / 3 m, 2 4 / 3 *) . (131) 

Interestingly, it turns out that the integral over z in Eq. (I126P can be performed explicitly, yielding 

P(s, «,) ~ ± I 2« (. + „*) (4, + „*) Ai (^^ j - 2'»/» (., + u?) Ai' (^^if N ' 

which allows to obtain straightforwardly the large \w\ behavior of P(s,w), valid for s — > oo 

logP(s, to) ~ ~\w\* ~l\w\s + O(log \w\) . (133) 

Note that the lead ing term matches perfectly with the behavior obtained in the large deviation regime, 
when M -> ^/2N 1(73)1. with w -+ 2 7 / 3 i;. 

^.5.3 Marginal probability density function of the position of the maximum and its large w behavior 

We now investigate the marginal pdf P(w) of the position of the maximum, which corresponds to the 
distribution of the end-point X/T 2 / 3 of the directed polymer (Fig. [T]). It is obtained from the jpdf 
P(s,w) by integration over s: 

P H = ^7s/ dsT x {s) J f(x,w)f(x,-w)dx, (134) 

which is obviously symmetric P(w) = P(—w), as it should be. It is normalized, thanks to the initial 
formula we started with ((2"T)) , though we could not check it directly from the above formula (|134[) . 
The distribution of the endpoint X of the directed polymer of length T has been widely studied in the 
physics literature. It was argued from scaling argument as well as from the study of simplified model, 
the so called "toy model" [531155] . that \ogP(X) oc ~(X/T 2 / 3 ) 3 for X > T 2 / 3 . The expression above 
(|134p yields an exact expression for the distribution of XjT 2 ! 3 ((5) [TTj) from which we can extract the 
asymptotic behavior. 

To analyze the large w behavior of (|134p the explicit formula obtained above (|125[) is very useful. 
Indeed, in the limit w — > +oo, the two quantities f(x,w) and f(x,—w) behave quite differently. For 
f(x,w), the integral over £ is dominated by the region of small £, due to the factor e~ w ^ while 
f(x, —w) is instead dominated by the region of large £, due to the term e w< * . One can then obtain the 
leading behaviors of f(x,w) and f(x,—w) by estimating the small and large £ behavior, respectively, 
of <p2(C,x). The details of this analysis are left in Appendix (A) It yields the asymptotic behavior 

logP(to) = -j^w 3 + o(w 3 ) . (135) 

Finally, from Eq. (|131[) one obtains the asymptotic behavior of the pdf P(t) of the position of the 
maximum of the Airy2 process minus a parabola (0J as 

\ogP(t) = ~t 3 + o(t 3 ), (136) 

which is in agreement with the result announced, without any detail, in Ref. [T7], where the amplitude 
of the cubic term was however not computed. Note also that this exponent — |i 3 was very recently 
given as an upper bound on the decay of the tail of P(t) in Ref. |74| . 
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5 Conclusion 

To summarize, we have studied the joint probability density function Pn(M,tm) of the maximum M 
and its position tm (see Fig. \2§ for N non-intersecting excursions in the large N limit. Once M and 
tm are properly shifted and rescaled we have shown that the jpdf converges, when A^ — ¥ oo, to a jpdf 
P(s, w) (jllOp which we have computed explicitly (|111[) . This jpdf P(s, w) yields, up to a rescaling f|131[) . 
the jpdf of the maximum and its position of the Airy2 process minus a parabola ^. The expression 
obtained here is actually different from the one (J127I) recently obtained in Ref. [17J. In particular we 
show here that this jpdf can be expressed in terms of the psi-function for the Hastings-McLeod solution 
to the Painleve II equation. After this paper has been submitted, the equivalence of the two formulas 
was shown in Ref. [75] . 

From P(s,w) one obtains an exact expression for the marginal pdf P(w), from which we have 
obtained the asymptotic behavior log P(u>) ~ —w 3 /12, establishing on firmer grounds this cubic be- 
havior which had been proposed some time ago in the physics literature, based on scaling argument 
PQ, the analysis of the "toy" model [53] and numerical simulations [571158) . Very recently, this pdf was 
also measured experimentally, indicating also such a cubic behavior in the tail of the pdf, while a 
more precise comparison between our exact result and the experimental data would certainly be very 
interesting. 
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A Asymptotic behavior of the marginal probability density function P(w) 

In this appendix, we study the asymptotic behavior of the (marginal) pdf P(w) of the position of the maximum. 
The starting point of our analysis is the formula given in the text in Eq. (11341) . which we recall here 

2 /*oo /'oo 

P M = ^2073 / dsTi(s) f(x,w)f{x,-w)dx . (137) 

We first obtain the leading behaviors of f(x,w) and f(x,—w) (|125|) by estimating the small and large £ 
behavior, respectively, of 4>a (£,#). 

One can obtain the small £ behavior of <&2{(,,x) by analyzing the coupled equations for $i((,,x),$2(C,,x) 
in (I112II113|) and using that #i(— C,,x) = $i((,,x) while #a(— C,,x) — — ^a(C^)- One obtains 

$i(C,x) = Ce- I ~ 9< - v)dv + 0(C i ) , (138) 

4> 2 (<;,x) = -Ce f Z° q(v)dy [C f X dze- 2f ^ q(y)dv -D] + 0(C 3 ), (139) 

J — OO 

where at this stage the constants C and D remains undetermined. The large w > behavior of f(x,w) is then 
given by 

nl3/3 roo 2 rjl3/3 rx 

f(x, W ) = -=-=- / ^ 2 ^,x)e' W< - d( ~ -f—j-ef* dy q{ y)dyp I ^ g - 2 /, q (y)dy _ ^ _ (14Q) 

The constants C and D can then be obtained by extracting the large x behavi or of the above expression (|140|) 
and matching it with the large w expansion of the expression given in E g. (I104|) . assuming that the limits 
w — > co and x — > oo do commute. When x — > co the integral over z in Eq. (|140[) is diverging and it is easy to 
see that one has 

C ^ dze- 2f ~ q(y)dv ~Cx, x>l. (141) 

■J — oo 

Therefore, to leading order in x, for x — > oo, one obtains from Eqs. (|140|) and (|141|l 

2 7/3 Cx 
/(-'») ~-5-i7T, (142) 
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while the expression in Eq. (|104|l yields 

/(*.«o = JS^ + °(-" 5/2 ). ( 143 ) 

which, by comparison, yields immediately 

C = 0F . (144) 

The computation of D is slightly more involved but by comparing Eq. (|140l) in the large x limit one obtains 



D = C 



/0 roo 

dze~ 2 / ~ ^'^ + / dz(l - e" 2 J ™ 9(!/)d!/ ) 
-oo JO 



(145) 



Given that q(y) > 0, for all y real, one sees immediately on that expression (|145[1 that D > 0. 

To compute the asympto tic b ehavior of fix, —w) for w — > oo we notice that, in this case, the behavior 
of the integral over C, in Eq. (|125[) is instead dominated by the regio n (" — » oo where, to leading order in w, 
<?2(C £) can thus be replaced by its asymptotic behavior given in Eq. (|118p . Performing the integral over £ one 
thus arrives at the expression obtained previously in Eq. (|104[) . with the substitution w —¥ —w from which one 
has to extract carefully the large w behavior. One obtains 

o7/3 3 

/(a,, - w ) „ ___ w VVn— r . (146) 

These asymptotic behaviors in Eqs. (|140l I146|) suggest to separate the integral over s in Eq. (| 13Tf) into two 
parts, 

P(w) = P+{w) + P-(w) 

2 roc roo 

P +H = 52073/ dsTx{s)J f(x,w)f(x,-w)dx, (147) 

2 ,0 /-oo 

P -H = 52o73/_ dsT4s)J f(x,w)f(x,-w)dx. (148) 

(149) 

Let us first analyze P+(w) in Eq. (I147JI . There, because of the exponential term e~ w x ' 2 coming from (|146[1 one 
can perform, in the integral over x, the change of variable z = vox to obtain 

c 1 3 « 
f( x , w )f( x ,- w )dx~-^e-T2™ ~- , (150) 

where c can be read from Eqs. (J14Q|[T46)} : 

q9/2 />oo 

c=^— / dz(l- e - 2 ^ 9(!/)d!/ ), (151) 

I"' 1 Jo 

where we have used the explicit expression of C in Eq. (|145p . Finally, from Eq. (|150j> one obtains the large to 
behavior of P+(w), to leading order as 

Pi 1 o 9 — 7/6 r+00 

P + ( w ) ~ ±_ e -T5- c , = - J"i(0) / d«(l-e- 2/ - «(»)<«»). (i 5 2) 

Let us now analyse P-(w) in Eq. (|148|) . whose analysis turns out to be more complicated. W e first n otice 
that J-i(s) is bounded for s < and in addition, its asymptotic behavior for s — > —00 is given by [761177] 

o-H/48 , _|,|3 |a| 3/2 

■ p i( s ) = ^T7Tr eC( " 1)/2e ^^"^^( 1 +°(i s r 3/2 )). ( 153 ) 

so that there exists a constant K such that 

Ji(s) <Ke~ s2 , Vs<0. (154) 

Therefore one has 

/o 2 r 00 

dse"" / f{x,w)f{x,~w)dx . (155) 

-OO «/ S 
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For large w, because of the exponential term e~ wx ' 2 in Eq. (|146|) . the integral over s in Eq. (|155|l is dominated 
by the region of large negative s. On the other hand, for large negative s, the integral over x is also dominated 
by the region where x is large and negative. Now using the behavior of q(s) (see e.g. Ref. |76] ) 

/ — s 
?(«) ~ y — , s-¥ -oo , (156) 



one obtains that, for large negative s, one has 

f(x,w)f{x,-w)dx~ — e"H*'" + " J 3 ^ , (157) 



/' 

J S 



d _ i m 3 + v^| a j 3 / 2 . 



where d > is a constant which can be read fr om E qs. 1)14011146)1 . Using this last estimate (|157|l . one immediately 
sees that the large w behavior J_ (w) in Eq. (|155p is given by 

log I - (w) = --^w 3 + o{w 3 ) . (158) 



Finally, using the estimates in Eq. (|152|l and (|158[) together with the inequality in Eq. (|155[) one obtains 



as given in the text 1)135 



log P(™) = --^-w 3 + o(w 3 ) , (159) 
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